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Relativistic Hartree-Fock method together with many-body perturbation theory and configuration 
interaction techniques are used to calculate relativistic energy shifts for frequencies of the strong 
electric dipole transitions of C III, C IV, Na I, Mg I, Mg II, Al II, AI III, Si IV, Ca II and Zn II. These 
transitions are used for search of the variation of the fine structure constant in quasar absorption 
spectra. The results are in good agreement with previous calculations. The analysis of Breit 
contributions is also presented. 
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I. INTRODUCTION 

Search for variation of fundamental constants is cur- 
rently an extensive area of research which spans the whole 
lifetime of the Universe from Big Bang nuclear synthe- 
sis to present day atomic-clock experiments (see, e.g. a 
review by Uzan 1]). This search is motivated by theo- 
ries unifying gravity with other interactions as well as by 
many cosmological models. In both cases a possibility 
for fundamental constants to very in space and/or time 
is present. 

Strong evidence that the fine structure constant, a 
might be smaller about ten billion years ago was found in 
the quasar absorption spectra [3, H, 0, |^,H, 0| ■ This re- 
sult was obtained from the analysis of the data from Keck 
telescope in Hawaii by the group of researchers based at 
University of New South Wales in Australia. However, 
the analysis of the data from VLT telescope in Chile per- 
formed by different groups ^, i^] gives null result. There 
is an outgoing debate in the literature about possible rea- 
sons for this disagreement. 

All the analysis in Refs. 0, i, i, i, i, 0, i, i] was 
performed with the use of the so-called many-multiplet 
(MM) method which was first suggested in Refs. 0, ■ 
This method uses frequencies of strong atomic electric 
dipole transitions for the analysis. Its sensitivity to vari- 
ation of the fine structure constant is more than an order 
of magnitude better than the analysis of the fine struc- 
ture intervals which was used before [111, [Hi [ll] • This 
dramatic gain in sensitivity comes with some complica- 
tions. The method relies on atomic calculations to reveal 
dependence of atomic frequencies on the fine structure 
constant. All calculations used in the analysis so far were 
performed within a single group of researches bases at the 
University of New South Wales [H, [H, [H, [13, [S [H ■ 



Due to importance of detecting of any variation of fun- 
damental constants and controversy of the results it is 
important that atomic calculations be also verified by 
independent calculations. 

A positive development in this direction is the recent 
independent calculation of the relativistic energy shifts 
in Fe II ion . Fe II is a single most important element 
for the analysis of quasar absorption spectra. It has lines 
which move in opposite direction if alpha varies, and the 
value of this shift is relatively large. In principle, sub- 
ject to sufficient statistics, it alone can serve as a probe 
of variation of the fine structure constant in quasar ab- 
sorption spectra j20j. However, calculations for Fe II are 
difficult due to large number of valence electrons. 

In present work we perform further verification of the 
relativistic energy shifts in atoms of astrophysical interest 
by considering atoms and ions with one and two valence 
electrons. In the case of single- valence electron atoms we 
use Dirac-Hartree-Fock method (DHF) and relativistic 
many-body perturbation theory (RMBPT). For double- 
valence electron atoms we use DHF, RMBPT and con- 
figuration interaction (CI) technique. The results on this 
stage of calculations are in very good agreement with 
previous calculations. 

We also perform the analysis of role of Breit interac- 
tion. It turns out that Breit contribution to the energies 
and relativistic energy shifts are small. However, Breit 
interaction gives significant contribution to the fine struc- 
ture intervals brining them to excellent agreement with 
experiment. Therefore, inclusion of Breit contributions is 
important for the analysis of the accuracy of calculations. 



II. METHOD 



A. q-factor calculations 
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The difference between frequencies in QSO spectra and 
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in the laboratory after taking into account the Doppler 
shift depends on values of fine-structure constant a. For 
small changes of a, the transition frequency changes lin- 
early with (a/ao)^ and can be presented in a form 



uj{x) = ujo + qx, 



(1) 



where luq is the laboratory value of the frequency and 
X = (a/ao)^ — 1, q is the sensitivity coefficient to be 
found from atomic calculations. Note that 
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dw 
dx 



at X — 0. 



(2) 



To obtain q-factors theoretically, it is necessary to calcu- 
late energies for at least two different values of the fine 
structure constant. Symmetric formula for the deriva- 
tive is more accurate, and better precisions of q can be 
reached by calculating energies at two values of alpha 
symmetrically displaced from the standard value. Due 
to numerical accuracy issues and non-linear dependence 
in a large range, the displacement should be not very 
small and not very large. We find that ai — 1/134.0 and 
a2 — 1/140.0 satisfy both conditions. The q values are 
calculated using the following equation 



q^[E{a,)^E{a2)]al/{al~al), (3) 
The energies and q-factors are both in inverse cm units. 

B. DHF and 2nd-order RMBPT 

For g-factors of univalent atoms and ions, we calculate 
energies using the 2nd-order RMBPT formalism built on 
expansion in the frozen-core 

^(A'-i) DHF basis. By com- 
parison of q-factors from the 2nd-order and DHF calcula- 
tions, we find that the precision of second-order theory is 
expected to be sufficient. The 2nd-order RMBPT formal- 
ism is described in Ref. [2l[ . In the DHF basis first-order 
correction is already included, and the number of second- 
order corrections is reduced. The summation over excited 
states in the second-order expressions is carried out by 
using a finite compact B-spline DHF basis, with cavity 
sizes chosen to minimize influence of boundary conditions 
on the valence energies of interest. Angular orbital mo- 
mentum of excited states is limited to 5, without much 
reduction in the precision. Other parameters are chosen 
to minimize numerical errors. 

The second-order RMBPT gives much more accurate 
energies than the DHF theory as can be seen from our 
calculations presented in Table [D We also compare theo- 
retical and experimental fine-structure splitting between 
P3/2 and pi/2 states. The agreement is very good in the 
second-order of RMBPT. However, it is further signifi- 
cantly improved when the Breit corrections discussed in 
next section are added. This is another indication of high 
accuracy of the calculations. 



TABLE I: DHF and 2nd-order energies and fine-structure 
splittings of univalent atoms/ions 



Element State 



DHF 



Energy 
E2" 



Fine Structure 
Expt' E2" E2+Br" Expt*- 



C IV 



Na I 



Mg II 



Al HI 



Si IV 



Ca II 



Zn II 



2pi/2 

2P3/2 
3pi/2 
3p3/2 
4pi/2 
4p3/2 

3pi/2 

3p3/2 
4pi/2 
4p3/2 
3pi/2 
3P3/2 
4pi/2 
4p3/2 
3pi/2 
3p3/2 
4pi/2 
4p3/2 
4pi/2 
4p3/2 
4pi/2 
4p3/2 



65201 


64548 


64484 










64680 


64592 


132 


107 


108 


15921 


16812 


16956 








15937 


16831 


16973 


18 


17 


17 


28904 


30068 


30267 








28909 


30074 


30273 


5 


5 


Q 


34530 


35603 


35669 








34620 


35700 


35761 


97 


91 


92 


78574 


80463 










78605 


80496 


80650 


32 


30 


31 


52709 


53672 


53683 








52944 






947 




234 


141252 


143538 


143633 








141334 


143623 


143714 


85 


80 


80 


70540 


71309 


71288 








71008 


71794 


71749 


486 


461 


461 


215704 


218226 


218267 








215870 


218397 


218429 


171 


162 


162 


23403 


25490 


25192 








23603 


25722 


25414 


232 


225 


223 


44610 


48548 


48481 








45347 


49429 


49355 


881 


866 


874 



''DHF-(-2nd-order 
''NIST, Ref. jsJl 
<=E2-|-Breit (see Table lUll 



C. Breit corrections 

Relativistic energy shift which was considered above 
is due to the difference between Dirac and Schrodinger 
equations. This difference leads to a correction to the 
energy proportional to in the leading order. There- 
fore, for small change of a this correction coincides with 
the definition of the q-coefficient (see formula (HJ . How- 
ever, there is also Breit relativistic correction to the 
inter-electron interaction 



22l |. This correction is also 



proportional to and therefore contributes to the q- 
coefficients. It is important to check the values of these 
corrections to have reliable results. We include Breit in- 
teraction using the technique developed in our previous 
works [23, [11]. 

We use the following form of the Breit operator (atomic 
units) 



Qi ■ 0:2 + {&! ■ n){a2 ■ n) 
2^ ■ 



(4) 



Here r = fir, r is distance between electrons and cti is 
the a-matrix of the corresponding electron. This is a 
low frequency limit of the relativistic correction to the 
Coulomb interaction between electrons. It contains mag- 
netic interaction and retardation. 

Similar to Coulomb interaction, Breit interaction cre- 
ates a potential which is to be added to Hartree-Fock 
potential. In the case of closed-shell atoms direct term 
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in Breit potential vanishes and only exchange term re- 
mains. This is the case for single- valence-electron atoms 
considered in present work since we use the T^^"^ ap- 
proximation. 

Self-consistent calculations arc performed for a closed- 
shell core in a potential which is a sum of Coulomb and 
Breit terms 

V^V^ + V^. (5) 

States of valence electrons are calculated in the same 
potential ([5|). In this approach Breit interaction be- 
tween electrons receives exactly the same treatment as 
the Coulomb one. It is first included as interaction be- 
tween core electrons and then as an interaction between 
valence and core electrons. Therefore, an important ef- 
fect of core relaxation is included. A less important effect 
of Breit interaction on inter-electron correlations is not 
included in present work. This is justified by small value 
of the corrections. 

Note that non-perturbative treatment of Breit interac- 
tion leads to inclusion of higher-order in Breit operator 
terms, terms proportional to (H^)^ ,{H^)^ , etc. Inclu- 
sion of these terms cannot be justified and in principle 
they can be easily eliminated by a rescaling procedure in 
which Breit operator is suppressed in the calculations by 
a scaling parameter A and then the answer is interpolated 
to A = 1. It turns out, however, that as a rule A = 1 is 
already in linear regime. 

The values of the corrections are found by running pro- 
grams with and without corresponding extra terms in 
the potential. The calculated values are presented in Ta- 
ble ini In Table U we also included Breit corrections to 
the fine structure splitting. 

Calculations show that although first-order valence 
Breit correction is a dominant contribution among Breit 
corrections to the energy of valence electrons, transition 
energy has a substantial cancelation for this correction, 
and Breit core-relaxation contribution becomes compara- 
ble with the valence Breit contribution for the transition 
energies and hence for q- values. The inclusion of core- 
relaxation effect significantly changed the value of Breit 
contribution. 



D. BO+CI method 

To calculate alpha variation coefficients for diva- 
lent atoms and ions, we will use the Brueckner-orbital 
(BO)-fCI method, introduced in Ref.jl^, and modified 
in Ref.|2(|] to include first-order Breit corrections. This 
method is essentially the combination of CI, to treat 
strong valence-valence interactions, and MBPT, to treat 
important valence-core interactions. It is similar to the 
method discussed in Refs.[27l. l28l. [29|. 

The BO-CI method described in detail in Ref.[25t 
is based on the effective Hamiltonian formalism which 
leads to the problem of diagonalization of the Hamil- 
tonian matrix built on the two-electron configuration 



TABLE II: Calculations of g-factors for univalent atoms/ions 



Atom/Ion 


State 


DHF+2nd Breit 


Total 


Other" 


C IV 


3pl/2 


102 


13 


115 


104(20) 




3p3/2 


233 


-12 


221 


232(20) 


Na I 


3pl/2 


45 


-1 


44 


45(4) 




3p3/2 


63 


-2 


61 


63(4) 




4p3/2 


59 


-2 


57 


59(4) 




4pl/2 


53 


-2 


51 


53(4) 


Mg II 


3pl/2 


119 


1 


120 


120(10) 




3p3/2 


216 


-5 


211 


211(10) 




4pl/2 


167 


-6 


161 






4p3/2 


200 


-8 


192 




Al III 


3pl/2 


218 


5 


223 


216(14) 




3p3/2 


466 


-9 


457 


464(30) 




4pl/2 


349 


-12 


337 






4p3/2 


434 


-17 


417 




Si IV 


3pl/2 


347 


13 


360 


346 




3p3/2 


835 


-12 


823 


862 




4pl/2 


617 


-20 


597 






4p3/2 


789 


-29 


760 




Ca II 


4pl/2 


219 


3 


222 


224 




4p3/2 


454 


-4 


450 


452 


Zn II 


4pl/2 


1590 


-5 


1585 


1584(25) 




4p3/2 


2508 


-20 


2488 2479(25) 



'Ref. 113 



TABLE III: Calculations of g-factors for divalent atoms/ions 

Atom/Ion State CI-HVIBPT Breit Total Other" 
"cm 2s2p 163 165 

Mg I 3s3p 93 -7 85 86(10) 

3s4p 89 -8 80 87 
Al II 3s3p 270 270(30) 

"Ref. |3Q1 

state functions. Beyond the frozen-core Hamiltonian the 
first-order electron-electron interaction Hamiltonian and 
second-order correction which consists of the two-particle 
screening correction and the one-particle self-energy cor- 
rection are included. In the BO-CI method, the basis 
functions are chosen as BO and include second-order self- 
energy corrections together with DHF potential. The 
residual two-particle Hamiltonian matrix, that includes 
first-order valence-valence interaction and second-order 
Coulomb screening interaction, is evaluated in the BO 
basis and diagonalized to obtain state energies and CI 
wave functions. 



III. RESULTS 

Results of calculations of g-factors for univalent and 
divalent atoms and ions are presented in Tables |TT] and 
lllli respectively. Accurate agreement is achieved between 
our second-order values and g-factors previously reported 
and compiled in Ref. (soj . Univalent atoms and ions are 
calculated from second-order RMBPT energies with the 
method described in the previous section. Because the 
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contribution from the second order turned out to be rela- 
tively small, we expect that second-order results will give 
quite reliable values. Because previously only dominant 
relativistic effects were included within DHF formalism, 
we also added Breit corrections to energies to investigate 
the effects beyond Dirac-Fock approximation. 

For divalent atoms and ions calculations are performed 
with the BO+CI code, which is described in the previous 
section. All first-order Breit corrections introduced into 
CI-KMPBT in Ref.^ are also included. 

IV. CONCLUSIONS 

We have calculated ^-factors for mono- and divalent 
atoms and ions of interest for the extraction of fine- 



structure variation from quasar spectra. Our results 
agree with good accuracy with previous calculations and 
provide necessary independent verification. In particular, 
more difficult for theory divalent atoms and ions are cal- 
culated with a new method, BO+CI. Breit corrections, 
ignored previously, have been also evaluated. Although 
found in this work to be small, potentially they consti- 
tute a dominant class of relativistic corrections beyond 
the Dirac-Fock formalism. 
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